Off-Diagonal Long-Range Order in Bose Liquids: Irrotational Flow and Quantization 

of Circulation 
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On the basis of gauge invariance, it is proven in an elementary and straightforward manner, but 
without invoking any ad hoc assumption, that the existence of off-diagonal long-range order in 
one-particle reduced density matrix in Bose liquids implies both the irrotational flow in a simply 
connected region and the quantization of circulation in a multiply connected region, the two fun- 
damental properties of a Bose superfluid. The origin for both is the phase coherence of condensate 
wave-functions. Some relevant issues are also addressed. 
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The discovery of superfluidity in liquid ^He leads to 
London's proposal that the transition from He I phase to 
He H phase is an analog of Bose-Einstein condensation 
(BEG), namely, the superfluidity in He H is characterized 
by the macroscopic occupation of a single quantum state 

0, despite that the liquid He II is actually a strongly 
interacting Bose system. On the basis of London's idea, 
Tisza 1^ , vi^hile Landau independently |§] , developed phe- 
nomenologically a two-fluid model, in which the super- 
fluid He II is regarded as being composed of superfluid 
and normal fluid components. The superfluid compo- 
nent is the single condensed quantum state, carrying no 
entropy, and thereby resulting in the assumption that 
the superfluid flow is irrotational. Meanwhile, Onsager 
Q suggested that quantized vortex lines, or the quanti- 
zation of vorticity according to Feynman independently 
1^, must exist in He II, which was experimentally veri- 
fied a few years later ||^. Since then, the fact that the 
irrotational flow and the quantization of circulation (vor- 
tices) are two fundamental characters of superfluid He II 
is well established. Recently, the existence of vortices in a 
coherent collection of ultracold atoms has also been pro- 
posed and experimentally observed by spinning up the 
condensate Q]. On the other hand, Penrose and Onsager 
1^ generalized the mathematical description of BEG to 
interacting Bose systems, and manifested that BEG ap- 
pears as the off-diagonal long-range order (ODLRO) in 
one-particle reduced density matrix exists, and liquid he- 
lium II was shown to possess such an order. Later, the 
concept was successfully extended to superconductivity 
in interacting fermion systems by Yang It is now 
well accepted that both superconductivity and superflu- 
idity are quantum phases characterized by the existence 
of ODLRO. Nevertheless, what is the relationship be- 
tween the fundamental properties of the condensates and 
the existence of ODLRO? The answer to the question for 
superconductors was presented a few years ago ]lO| , pd| , 

1. e., the hypothesis of ODLRO in the two-particle re- 
duced density matrix implies both Meissner effect and 



flux quantization, while the proof for superfluids is still 
absent. 

In this paper, by noting the analog between supercon- 
ductors and superfluids, we shall extend the proofs for 
superconductors presented in Refs. |l^,|lll to interact- 
ing Bose systems, and substantiate that the existence of 
ODLRO in one-particle reduced density matrix in a Bose 
liquid implies both irrotational flow and quantization of 
circulation. The origin for the two fundamental proper- 
ties of a Bose condensate comes from the phase coherence 
of condensate wave- functions. Our proof does not invoke 
any ad hoc assumption. 

Let us consider a bucket of Bose liquid composed of N 
homogeneous interacting spinless particles, rotated with 
a constant angular velocity ft. The Hamiltonian of the 
system is 



H 



mVsir,)r + lY.^ir,i), (1) 



where pj 



—ih\7j is the momentum of jth particle, 
V{rji) denotes the interactions between particles, and 
Vs {vj ) represents the drift velocity, given by 



Vs(r) = r2 X r + V6'(r) 



(2) 



in which the first term comes from the Coriolis force, 
and the second term, to be identified below, originates 
from the gauge invariance. Eq. comes out from 

the basic relation n =iV x v^, as is in a classical 
case. The Schrodinger equation reads H^pn{ri, • • • , tat) = 
Enipn(ji, ■ ■ ■ , rw), where the orthonormal eigenfunctions 
tpnS are single- valued and symmetric. In analogy to the 
spirit of Refs. |l^Jll[], we first make a gauge transfor- 
mation. Under a small spatial displacement r — > r — 1, 
we have Vs(r)— > Vs(r — 1) = Vs(r) + V[l • {ft x r)-|-6'(r — 
1) — 9{r)], and V{rji) — ^ V{rji). Then, the Schrodinger 
equation becomes 



, Tat) = Enlpniri, ■ ■ • , Tat) 



(3) 



where 



1 



V-nCri, • • • , ta,) = ^^'"^■''V«(i-i - 1, • • • , rjv - 1) 

(4) 

is also the eigenfunctions of H, and the phase factor 
^(r, 1) is defined by 



C(r, 1) = 1 ■ (O X r)+e'(r - 1) - 0(r). 



(5) 



Therefore, we have another set of orthonornial eigenfunc- 
tions, -iAn's, for the system. 

To gain insight into the underlying physics behind 
ODLRO in Bose liquids, we ought first to look at the in- 
trinsic properties of the one-particle reduced density ma- 
trix pi(r,r'). In accordance with the definition, pi(r, r') 
can be written down as 



Pi(r,r') 



dY2 ■ ■ -dYM 1 

(iV-1)! Z ^ 
'0„(r,r2, • • • ,rAr)?/'*(r',r2, • • ■ ,yn), (6) 



where Z — partition function. 

pi (r, r') is normalized as Trpi — J drpi (r, r) = N, and is 
single-valued. Now, we utilize the energy eigenfunctions 
given by Eq. to represent pi{r, r') anew, and get the 
relation 



pi(r, r') = e'xK(r.i)-«(rM)]^^(j. _ r' - 1). 



(7) 



Equation illustrates that under a spatial displacement 
the one-particle reduced density matrix gains a position 
and displacement dependent phase factor. It is an in- 
trinsic property of pi(r,r'). On the other hand, we can 
express pi(r,r') in the spectral form as 



pi(r,r') 



:A,$,(r)ci>:(r'), 



(8) 



where $1, is the eigenfunctions of pi, with eigenvalues A^. 
According to Penrose and Onsager the existence of 
ODLRO in pi implies the spectral resolution such that 



pi(r,r') = Ao$o(r)$S(r')+p'i(r,r'), 



(9) 



where Aq, on the order of 0{N), is the largest eigenvalue 
of pi, $0(1") is the corresponding eigenfunction (the so- 
called condensate wave- funct ion ) , and p[{r,r') denotes 
other terms being on the order of 0(1). In the thermo- 
dynamic limit, p[{r,r') as |r — r'| — > 00. Com- 
paring Eqs. M) and (Bh we may obtain under the limit 



r — r 



00 the following expression 



$o(r)$5(r') = e^f- -«('-M)l$o(r - l)%{r' - 1). (10) 

With a careful analysis, one may find that this equation 
suggests 



$o(r) =e^«(')e'^«(""'')$o(r-l), 



(11) 



where a displacement dependent, but position indepen- 
dent phase factor, C,{1) with C,{0) = 0, may possibly 
emerge when extracting the information of individual 
function from the products of Eq. (p^. As a result, 
we observe that the condensate wave- function, <I>o(r), is 
mutually related by a phase factor at different positions. 
It is this property, being of fundamental importance, that 
enables us to prove the irrotational flow and the quanti- 
zation of circulation in Bose superfluids. 

In a simply connected region inside the superfluid, 
according to the rule governed by Eq. (^fj) we per- 
form two opposite displacements, first a followed by b, 
and then b followed by a. Through comparing the 
relations satisfied by <&o(r) in the two operations, we 
discover that the phase factors satisfy the equation: 
F(r, a, b) = F(r, b, a), with the function i^(r,a, b) — 
exp{i^[^(r,a) -I- ^(r — a, b)]}. By solving this equation, 
we have exp[z^r2 • (b x a)] = 1, yielding 



n 



(bxa)=n— , 



(n = integer). 



(12) 



It is seen that the right-hand side of Eq.(|l2|) is dis- 
cretized, while the left-hand side can be continuously 
changed as the two operations can be arbitrarily choosen. 
The consistency condition requires 



n = o. 



(13) 



This proves the property of the irrotational flow of the 
Bose condensate. However, we should stress that in terms 
of Eq. (||) the flow velocity, i.e., the superfluid velocity, 
should be in principle nonvanishing, i.e., Vs(r) = V6'(r), 
satisfying V x Vs(r) = 0. 

In a multiply connected region inside the superfluid, 
there should be some singularities around which the 
fluid is rotating. These singularities are nothing but 
the vortex lines. In this case, we can repeat succes- 
sively for infinitesimal displacements along a closed path 
C which encloses at least one vortex line, and then pick 
up products of phase factors. Consequently, we have from 



Eq. (|l l[) the relation $o(r) = e ~ '''''"'$o(r), where use 
has been made of the single- valuedness of <I'o(r), be- 
ing a consequence of the single- valuedness of pi, and 
exp[iC(fp dl)] = exp[i^(0)] =1. Therefore, we get 



(f Vs ■ dl ~ n — , (n = integer). 
Jc 'm 



(14) 



This is the quantization of circulation, with the quan- 
tum of vorticity being h/m. If there is no any vortices 
encircled by the path loop C in the multiply connected 
system, then the circulation will be identically zero, lead- 
ing to again the property of irrotational flow everywhere 
away from the vortex cores. 

The irrotational flow and the quantization of circula- 
tion have been well understood for superfluid He II based 



2 



on the work of London, Onsager and Feynman quite long 
time before. Here we have offered another route, without 
ad hoc assumption, to understand the two fundamental 
characters of a Bose superfluid, and then established a 
deep connection between the basic properties of the con- 
densate and the existence of ODLRO in Bose liquids. 
It is manifested that the condensate, as a single quan- 
tum state possessing these two essential features, is in- 
deed characterized by the eigenstate of the one-particle 
reduced density matrix with the largest eigenvalue of the 
order of 0{N), as argued by Penrose and Onsager 
The preceding discussion shows that the phase coherence 
of condensate wave-functions is the origin for occurence 
of both properties. 

Next, let us clarify the physical meaning of the quan- 
tity 6{r). In general, the condensate wave- function is as- 
sumed to be a complex number, and is usually denoted 

by 



0o(r) =0o(r-l). 



(18) 



$o(r) = 0o(r)e 



ia{r) 



(15) 



where the amplitude (/>o(r) and the phase factor a{r) 
are real functions. Define the condensate density as 
po(r) = |*i'o(r)|^ = [00 (r)]^. Then, the superfluid current 
density at position r is given by Jo(r) = --^{^qV^o — 
<i>oV<i>Q) = ■^[(j)o{r)]^ya{r). Since the current den- 
sity and the superfluid velocity is related by Jo(r) = 
po(r)vs(r), combining the irrotational property of the 
condensate we have 



e{r) = -a(r). 



(16) 



As a result, we arrive at the conclusion that the superfluid 
velocity in the Bose liquids is proportional to the gradi- 
ent of the phase factor of the condensate wave-function, 
namely 



Vs(r) = — Va(r) 
m 



(17) 



a well-known result that was supposed for superfluid He 
II by London, Onsager and Feynman. However, our 
derivation shows that Eq. (17) holds true for any Bose 
system so long as the system possesses ODLRO in pi. 
Thus, the phase of the condensate is closely related to 
the gauge invariance. 

To explore the displacement dependent but position 
independent phase factor ('(I), let us compare Eqs. ( p^ ) 
and (|l|). By considering (H) and (|lf) we can get 
tan[2a(r)] = tan[2a(r + C(l) f 1 • (fJ x r)]. Solving 
this equation and noting (^(0) = 0, one may obtain 



C(l) = 2[a(r)-a(r- 



l)]--l.(nxr). 



This suggests that such a phase factor indeed exists if the 
condensate wave- function takes the form of (^5|) . Conse- 
quently, we have from Eqs. (^l|) and ( p^ ) 



It shows that the amplitude of condensate wave-function, 
under assumption of the form (p^, is displacement- 
invariant. This is convincing, because in superfluid He 
II, 00 (r) is supposed as the square root of uq , the macro- 
scopic occupation number of the condensate, which is 
certainly translation- invariant. Here we have proved that 
for any condensed Bose system the condensate density 
possesses a translational invariance. 

In the Bose condensate, the continuity equation be- 
comes 90o/9i + Vs • V0O + (l/2)0oV • Vs — 0. After inte- 
grating it, we have 



d_ 
di 



j dv<j)o{T,t) = dr0o(r,i)V-v,(r) 



drV- [v,(r)0o(r,t)]. 



In view of Gaussian theorem, the second term of the 
right-hand side of above equation should be a surface 
integral, and is vanishing in the present case. For a 
steady flow, V-Vs(r) — 0, the above equation tells us that 
/ (ir0o(r, t) is independent of time. If we define the spa- 
tial Fourier transform of 0o(r, t) as 0o(k, t), we know that 
limk^o 00(1^, t) is in fact independent of i. Consequently, 
the temporal Fourier transform of 0o (k, t) , defined by 
(/io(k, w), has a property: limk-to 0o(k, t^) ^ This 
implies that in the long wavelength limit there must be 
excitations whose frequency is zero in the Bose conden- 
sate. These massless excitations are nothing but Gold- 
stone modes. As one knows, in He II the Goldstone 
bosons are identified as second sound. Therefore, by 
means of Goldstonc's theorem one can conclude that in 
the Bose condensate the gauge symmetry is broken. 

It is worth noticing that the irrotational flow and the 
quantization of circulation in a Bose superfluid can also 
be derived starting directly from Eqs. (|l5|) and (p^, as 
was done in previous literature. However, if we take a 
careful look, we can observe that the prerequisite of such 
reasonings is that the condensate wave-function was sup- 
posed to take the form of Eq. ( p^. Contrary to the 
previous arguments in which Eq. (|15|) was only an as- 
sumption without proof, the present derivation docs not 
depend on what form of the condensate wave-function 
is taken, and only invokes the general principle of gauge 
invariance. In some sense, our result suggests implicitly 
that the condensate wave-function has the form of Eq. 

(H). 

To this end, since superfluids and superconductors 
share some common features, it would be interesting to 
compare the correspondence between characters and vari- 
ables of both, as summarized in the following table. 
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Supcrfluids 



Coriolis Force 
Angular velocity (17) 
Flow velocity (vs(r)) 
ODLRO in pi 
Irrotational flow 
(f2 = iV X = 0) 
Quantization of circulation 
(/v.-dl = nA) 

Broken gauge symmetry 
Phase coherence 
BEC state of bosons 



Superconductors 



Electromagnetic force 

Magnetic field (B) 
Vector potential (A(r)) 
ODLRO in p2 
Meissner effect 
(B = V X A = 0) 
Flux quantization 
(/A.dl = ni) 
Broken U(l) gauge symmetry 

Phase coherence 
Condensate of Cooper Pairs 



As is seen, the similarity between both systems is a con- 
sequence of ODLRO, or equivalently, of the broken gauge 
symmetry, which can be viewed as the profound underly- 
ing physics behind these two important systems in con- 
densed matter, as emphasized by many authors (see, e.g. 
Ref. for more references). 

We close this paper with a few remarks in order. 

(1) The present proof works for both nonintcracting 
and interacting Bose systems. Thus, for an ideal Bose 
gas, as there is an ODLRO in pi, the BEC state should 
also possess the property of irrotational flow and quan- 
tization of circulation. An ideal Bose gas is a superfluid, 
but it does not exhibit superfluidity, because it is unsta- 
ble against any motion of the system. It is also man- 
ifested by Landau's criterion for superfluidity, i.e. the 
critical velocity is zero. 

(2) For fragmented condensates, namely there are more 
than one eigenvalues of the one-particle reduced density 
matrix, pi, of the order 0{N), leading to two or more 
eigenstates of pi being macroscopically occupied [p^ , 
our proof seems not to be sufficient, though the forms 
like Eq. (11) for condensates are still a possible solu- 
tion. The fragmentation of condensates usually occurs 
in multi-component Bose systems, and cannot happen 
in single component systems, because it costs a macro- 
scopic extensive exchange energy. For non-fragmented 
systems our argument is closely related to the so-called 
phase locking, as discussed by Noziercs [T^ ]. 

(3) For the Bose systems with spin nonzero, the one- 
particle reduced density matrix, pi{r, a;r' , a') with a 
and a' spin indices, has a spectral resolution such that: 
pi(r,cr;r',cr') = Ao$o(r,cr)$5(r', cr') -I- p[{r, a-r' , a'), 
where p'l — > for |r — r'| oo. Obviously, the conden- 
sate wave- function, ^Q(r,a), is a spinor. In the absence 
of external magnetic fields, our reasonings remain simi- 
lar to those for spinless systems. However, if a magnetic 
field is applied to the system, owing to a so-called "local 
spin-gauge symmetry" [|lj] there is a vorticity induced by 
spatial variations of the magnetic field, and the relation 
like Eq. (H) no longer holds, suggesting that our proof 
appears to be invalid, which is a possible consequence of 
the fact that the time-reversal symmetry is violated in 
the presence of an external field. 



(4) For trapped alkali atoms like "^Li, '^^Na, ^^K, ^'^Rb 
etc, the signature of BEC is successfully observed at ex- 
tremely low temperatures (see, e.g. Ref. |l^] for a re- 
view). These trapped Bose gases have finite sizes and 
are inhomogeneous, which contain the number of atoms 
typically from a few thousands to several millions. There- 
fore, the thermodynamic limit is not exactly reached for 
trapped gases, and consequently, the concepts of broken 
gauge symmetry and ODLRO appear, strictly speaking, 
not to be applicable in these finite-sized systems. How- 
ever, though the phenomena occurred in trapped alkali 
atoms are not a thermodynamic phase transition in a 
rigorous sense, our preceding discussion might shed some 
useful light on it. For instance, the relation (|l^) satisfied 
by the condensate wave-function, <I>o(i"); which is nor- 
mally considered as the order parameter, could also be 
used to investigate the linear, not rotational, motion of 
the condensate to a proper level of approximation. Quite 
recently, the irrotational fiow and the quantization of vor- 
ticity in a BEC of ^^Rb atoms has been successfully ob- 
served at extremely low temperatures fl^ . 
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